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Telomeres are repetitive DNA structures that, together 
with the sheUerin and the CST complex, protect the ends 
of chromosomes. Telomere shortening is mitigated in stem 
and cancer cells through the de novo addition of telomeric 
repeats by telomerase. Telomere elongation requires the 
delivery of the telomerase complex to telomeres through a 
not yet fully understood mechanism. Factors promoting 
telomerase-telomere interaction are expected to directly 
bind telomeres and physically interact with the telomerase 
complex. In search for such a factor we carried out a 
SILAC-based DNA-protein interaction screen and identi- 
fied HMBOXl, hereafter referred to as homeobox telomere- 
binding protein 1 (HOTl). HOTl directly and specifically 
binds double-stranded telomere repeats, with the in vivo 
association correlating with binding to actively processed 
telomeres. Depletion and overexpression experiments 
classify HOTl as a positive regulator of telomere length. 
Furthermore, immunoprecipitation and cell fractionation 
analyses show that HOTl associates with the active telo- 
merase complex and promotes chromatin association of 
telomerase. Collectively, these findings suggest that HOTl 
supports telomerase-dependent telomere elongation. 
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Introduction 

Telomeres, the nucleoprotein structures at the ends of 
chromosomes, consist of 5'-TTAGGG-3' repeats bound by a 
dedicated set of proteins forming the shelterin complex (Palm 
and de Lange, 2008) . Three members of the six protein complex 
directly bind telomeric DNA and were the only direct telomere- 
specific binding proteins known so far: TRFl and TRF2 bind 
double-stranded DNA (dsDNA; Zhong et al, 1992; Bilaud et al, 
1997; Broccoli etal, 1997), whereas POTl binds single-stranded 
5'-TTAGGG-3' repeats (Baumann, 2001). This complex 
constitutively associates with telomeres and shields the ends 
of linear chromosomes from being recognized as a double- 
stranded break, thus protecting telomeres from end-to-end 
fusions (Palm and de Lange, 2008). While this solves the 
end-protection problem (de Lange, 2009), maintaining 
telomere integrity itself is of outstanding importance and 
major factors may have remained elusive. 

Telomere length homeostasis, a crucial process in stem cell 
biology, aging and cancer, depends on the equilibrium 
between telomere lengthening (in most cases due to telomer- 
ase activity) and shortening reactions (generally due to 
replication and controlled processing) (Jain and Cooper, 
2010). Telomerase is capable of adding telomeric repeats to 
chromosome ends de novo. The enzyme works as a 
ribonucleoprotein complex, which consists of a catalytic 
subunit with reverse-transcriptase activity (called TERT), 
and an RNA serving as the elongation matrix for telomeres 
(called TR or TERC) (Greider and Blackburn, 1989). While 
these two core elements are sufficient for telomerase activity 
in vitro, biochemical analyses have shown that in vivo 
telomerase resides in a large complex of about 1 MDa 
(Schnapp et al, 1998). Some additional components of this 
large multi-subunit holoenzyme complex have been 
identified. In particular, the core components of box H/ACA 
small nucleolar ribonucleoprotein particles (snoRNPs), DKCl 
(dyskerin), GARl, NHP2 and NOPIO, are part of the active 
telomerase complex, and are necessary for proper RNP 
assembly as well as for TERC stability (Mitchell et al, 1999; 
Wang and Meier, 2004). More recently, the ATPases RUVBLl 
and RUVBL2 have been identified as factors essential for 
holoenzyme assembly (Venteicher et al, 2008), and TCABl 
(WDR79/WRAP53), identified as a DKCl interaction partner, 
was shown to be required for proper localization of CAB box 
containing small Cajal body (CB) -specific RNPs (scaRNPs) to 
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CBs, including TERC, and is part of the active telomerase 
complex (Tycowski et al, 2009; Venteicher et al, 2009). 

The presence of major scaRNA processing and trafficking 
factors in the telomerase complex hints to an important 
aspect of telomerase ceh biology: the orchestrated maturation 
of telomerase and interaction with telomeres in the CB. 
Telomere maintenance by telomerase requires that both 
TERT and TERC are recruited from distinct subnuclear sites 
to telomeres during S phase (synthesis phase] (Tomlinson 
et al, 2008) . Like other scaRNAs, TERC contains a common 
CB-specific locahzation signal and accumulates in CBs (Jady 
et al, 2004; Zhu et al, 2004], where it is found together with 
TERT (Tomlinson et al, 2008]. In a cell cycle-dependent 
manner, telomerase-containing CBs are then recruited to 
telomeres, suggesting that CBs represent an enzymatic hub 
in which telomere elongation by telomerase takes place (Jady 
et al, 2006; Tomlinson et al, 2006; Cristofari et al, 2007]. This 
trafficking model is further supported by telomere elongation 
defects in the absence of TCABl or presence of dysfunctional 
TCABl, disrupting TERC accumulation in the CB (Venteicher 
et al, 2009; Zhong et al, 2011]. Nevertheless, so far it remains 
elusive how telomeres are recruited to CBs, how this selective 
interaction is regulated and what drives the conversion from 
telomeres in a closed state, in which telomerase has little or 
no access, to telomeres in an open, accessible state. 

Telomerase is usually limiting and, under physiological con- 
ditions, acts preferentially on short telomeres (Hemann et al, 
2001; Britt-Compton et al, 2009], due to a well-established 
negative feedback loop mediated in cis by TRFl and POTl, 
likely by hiding the 3'-overhang, which serves as a template for 
telomerase (Loayza and de Lange, 2003]. Indeed, diminished 
loading of POTl or expression of a dominant-negative version 
lacking DNA-binding activity leads to telomere elongation by 
telomerase, and in vitro experiments have shown that POTl is 
competing with telomerase for its substrate (Loayza and de 
Lange, 2003; Ye et al, 2004; Kelleher et al, 2005; Lei et al, 2005]. 
However, POTl also interacts with TPPl, and both proteins 
together promote telomerase activity in vitro (Latrick and Cech, 
2010]. Furthermore, TPPl has been shown to be required for 
the recruitment of telomerase to its substrate in vitro and to 
telomeric chromatin in vivo (Xin et al, 2007; Abreu et al, 2010; 
Tejera et al, 2010; Zaug et al, 2010; Zhong et al, 2012]. While 
TPPl has been proposed as a telomerase recruiter, it does not 
completely fit the definition, since it has initially been 
described as a negative regulator of telomere length in 
telomerase-positive cells (Liu et al, 2004; Ye et al, 2004], 
although these results may in part be attributed to secondary 
effects, such as the lack of POTl tethering to telomeres in the 
absence of TPPl. Based on these bivalent results, a ying-yang 
model for telomerase recruitment and activity control has been 
proposed for TPPl-POTl (Xin et al, 2007]. Indeed, recently, a 
specific patch of amino acids on the surface of TPPl, the TEL 
patch (TPPl glutamate (E]- and leucine (L]-rich patch], has 
been identified as crucial for the TPPl function in telomerase 
recruitment and regulation (Nandakumar et al, 2012]. Analysis 
of point mutations within the TEL patch demonstrated that the 
TEL patch is physically and functionally distant from the 
portion of TPPl engaged in end protection, separating these 
two functions (Nandakumar et al, 2012]. 

Regulating the amount of telomeric repeats added by 
telomerase also involves how long telomerase can act on a 
given chromosome end. Therefore, turning off telomerase 
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needs to be regulated in addition to the recruitment step and 
processivity control. In a recent study, Chen et al (2012] 
described the human CST (CTCl, STNl and TENl] complex 
as a terminator of telomerase activity (Chen et al, 2012]. CST 
competes with POTl -TPPl for telomeric DNA and is increas- 
ingly enriched on telomeric DNA during late S/G2 phase, 
correlating with the period in the cell cycle when telomerase 
action is terminated. In agreement with a suppression of 
telomerase action, depletion of any of the three CST complex 
members led to a steady increase in telomere length in a 
telomerase-dependent manner. The authors suggested that CST 
binds to telomerase-extended 3'-ends and thereby suppresses 
telomerase access and further elongation (Chen et al, 2012]. 
However, depletion of STNl has been reported not to affect 
telomere length in various telomerase-positive cellular contexts 
(Wang et al, 2012] and murine CTCl -null cells do not exhibit 
the reported telomere lengthening phenotype (Gu et al, 2012]. 
The proposed model is nonetheless appealing, as human CST 
subunits stimulate DNA polymerase a-primase (Casteel et al, 
2009] and therefore CST binding to telomerase-extended 
3'-ends could initiate a switch from telomere elongation to 
fill-in synthesis. This mechanism could provide an autonomous 
end point to telomerase action at single telomeres, ensuring 
that every telomere is extended by telomerase once and only 
once during every cell cycle (Chen et al, 2012]. 

We reasoned that the identification of novel telomere- 
binding proteins would be a first step to the identification 
of additional factors implicated in telomere biology. SILAC- 
based quantitative mass spectrometry (MS] has been adapted 
for DNA-protein interactions (Mittler et al, 2009] and has been 
used by us and others successfully to identify factors binding to 
particular functional DNA fragments (Markljung et al, 2009; 
Butter et al, 2010; Bartels et al, 2011; Butter et al, 2012] . For our 
purpose, we apphed this approach to telomeric DNA in order 
to screen for telomere repeat-binding proteins and identified 
the protein HOTl (HMBOXl; homeobox telomere-binding 
protein 1]. HOTl had previously been described as a putative 
transcriptional repressor based on reporter gene assays (Chen 
et al, 2006] and had been identified as a telomere-associated 
protein by the proteomics of isolated chromatin segments 
(PICh] approach (Dejardin and Kingston, 2009]. Here we 
demonstrate that HOTl directly binds to telomeric DNA, and 
characterize this binding in atomic detail by resolving a crystal 
structure of the HOTl homeobox domain in a cocrystal with 
telomeric DNA. In vivo, HOTl localizes to a subset of telomeres 
with a higher degree of HOTl -telomere association in cellular 
contexts of elevated telomere processing. In addition, we show 
that HOTl associates with the active telomerase complex and 
that HOTl is required for telomerase chromatin binding. These 
findings suggest that HOTl contributes to the association of 
telomerase with telomeres and telomere length maintenance in 
various cellular settings, and classify HOTl as the first direct 
telomere-binding protein that acts as a positive regulator of 
telomere length. 

Results 

Identification of HOTl as a direct telomere repeat- 
binding protein 

To identify telomere-binding proteins we used polymerized 
biotinylated double-stranded oligonucleotides of the telo- 
meric sequence (5'-TTAGGG-3'] and a scrambled control 
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sequence (5'-GTGAGT-3'], separately immobilized on para- 
magnetic streptavidin beads and incubated with heavy and 
hght SILAC-labelled nuclear extracts from HeLa cells, respec- 
tively. Specific binding of proteins is detected by incubation 
of cell lysates encoded by 'heavy' amino acids (^^N- and 
^^C-labelled Lys and Arg) with the bait sequence, while a 



control sequence is incubated with 'hght', nonlabelled amino 
acids. Specific binders display a differential SILAC ratio, 
whereas background binders have a 1:1 ratio. After mild 
washing, bead fractions were combined and captured pro- 
teins were analysed by quantitative, high-resolution MS (Cox 
and Mann, 2008) (Figure lA). 
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Figure 1 Detection of specific telomere-interacting proteins. (A) A schematic of the quantitative SILAC -based DNA interaction screen with 
DNA oligonucleotides containing either the telomeric repeat or a control sequence. Specific interaction partners are differentiated from 
background binders by a SILAC ratio other than 1:1. (B) MS spectra of representative peptides from the 'forward' pull-down experiment. The 
heavy peptide partners are easily detected (red dots), while the light partner is barely observable (blue dots) in the mass spectrum. (C) Two- 
dimensional interaction plot: known shelterin components cluster together with HOTl, demonstrating enrichment at the telomere sequence 
compared to the control sequence. (D) Summary of the MS data for HOTl and the core shelterin components from the SILAC-based 
DNA-protein interaction screens carried out with nuclear extracts derived from HeLa and murine ES cells. 
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We identified all the six core shelterin components with a 
SILAC ratio of about 10 or higher in the 'forward' and about 
0.1 in the 'reverse' experiment, in which we had switched the 
labels (Figure IB-D, Supplementary Figure SI and Supple- 
mentary Table SI). In contrast, none of the proteins known to 
interact with shelterin were identified with SILAC ratios 
sufficiently high to be consistent with telomere binding, 
demonstrating that this approach was very stringent and 
exclusively detected telomere repeat-binding proteins and 
their strong interaction partners (Figure IC). In addition to 
the shelterin components, we found the protein HOTl with a 
high SILAC ratio that clustered with those of the shelterin 
components (Figure IB-D, Supplementary Figure SI and 
Supplementary Table SI). This indicates that HOTl must 
either strongly associate with the shelterin complex or 
directly bind to the 5'-TTAGGG-3' repeats. 

To verify that the HOTl identification was not cancer-, 
cell- or species-specific, we repeated our telomere-binding 
assay with SILAC-labelled nuclear extracts derived from 
mouse embryonic stem cells (ES cells). Again, all compo- 
nents of the shelterin complex and HOTl were identified with 
SILAC ratios, indicating specific binding to the telomere 
repeats (Figure ID and Supplementary Table S2]. Here we 
also identified the two paralogues, POTla and POTlb, which 
result from a gene duplication of the Potl gene in the rodent 
lineage (Hockemeyer et aU 2006), underscoring the specificity 
of our assay for direct telomere-binding proteins. Hence, 
HOTl is a putative telomere repeat-binding protein 
conserved in mammalian cells. 

HOTl contains a homeobox domain (Chen et al, 2006], 
suggesting that it may bind DNA directly. To determine 
whether HOTl was detected in our assay due to direct 
binding to the 5'-TTAGGG-3' repeats, we performed DNA- 
binding assays with HOTl in vitro. Recombinant HOTl bound 
specifically to telomeric repeats, whereas no binding to the 
negative control repeat fragments (5'-GTGAGT-3') was 
detected (Figure 2A]. Exhibiting similar binding behaviour 
as TRFl, HOTl was not enriched on any of the subtelomeric 
variant repeats 5'-TCAGGG-3', 5'-TGAGGG-3' and 5'-TTG 
GGG-3', nor on the C. elegans telomere 5'-TTAGGC-3' repeat 
sequence (Wicky et aU 1996; Figure 2A). To test whether 
HOTl also associates with telomeres in vivo, we performed 
chromatin immunoprecipitation (ChIP] experiments with ex- 
tracts from HeLa cells using an antibody directed against 
HOTl. Similar to TRF2, HOTl IPs showed enrichment of 
telomeric DNA in comparison to two negative controls 
(anti-GFP antibody and IgG; Figure 2B]. Thus, HOTl is a 
direct and specific telomere repeat-binding protein. 



H0T1 recognizes telomeric DNA by means of its 
liomeodomain 

Driven by these findings and with the aim to fully understand 
the molecular interactions between HOTl and telomeric DNA, 
we crystallized the DNA-binding domain (DBD) of HOTl with 
telomeric DNA. In order to identify a construct suitable for 
crystallization, we initially tested six different HOTl fragments 
for their DNA-binding ability (Supplementary Figure S2]. 
The three longer constructs Q144-A345, L156-A345 and 
G233-A345 all bound to immobilized telomeric dsDNA baits, 
demonstrating that the homeodomain of HOTl is sufficient for 
recognizing telomeric DNA and that integrity of the predicted 
N-terminal POU-specific (POUs] domain (Chi et al, 2002] is not 
required. The three shorter constructs (P242-A345, P254-A345 
and R271-A345] were not able to bind to the bait DNA 
(Supplementary Figure S2] . For crystallization, we reconstituted 
and purified telomeric DNA complexes with all three binding 
constructs, but only one (G233-A345) yielded crystals when 
reconstituted with a duplex telomeric DNA (5'-cTGTTAGGGTTA 
GGGTTAG-3' and 3'-ACAATCCCAATCCCAATCt-5') similar to the 
one present in the crystal structures of the TRFl and TRF2 
homeodomains bound to telomeric DNA (Court et al, 2005] . 

The optimized crystals diffracted to 2.9 A resolution and 
we could solve the structure by molecular replacement using 
the NMR model of the human HOTl homeodomain (residues 
268-343, PDB entry 2CUF, unpublished data from RIKEN 
Structural Genomics/Proteomics Initiative). In the orthor- 
hombic crystals, the DNA forms an infinite double heHx via 
a C-T nonWats on- Crick base pairing of the single-base 
overhangs (not shown). Two copies of the HOTl DBD are 
bound to one duplex DNA that comprises two and a half 
5'-TTAGGG-3' repeats in a regular and undistorted B-form 
conformation (Supplementary Figure S3). As expected, and 
in accordance with the NMR model, the homeodomain of 
HOTl folds into a small structure of three consecutive helices, 
al (res. 276-288), a2 (res. 293-309) and a3 (res. 322-342), 
separated by a loop (al-a2) and a turn (a2-a3; Figures 2C 
and 3A, and Supplementary Figure S3). The N-terminal 
residues 233-266 and the two C-terminal residues were not 
defined by electron density and were, thus, not built. 

As reported for TRFl and TRF2, each copy of the HOTl 
homeodomain binds to the major groove of the DNA double 
helix around a 5'-TTAGGG-3' motif (Figure 2C and 
Supplementary Figure S3). Typical for homeobox domains, 
binding to DNA is mainly mediated by an N-terminal 
unstructured arm (267-276), the loop between al and a2, 
and the C-terminal a3 (Figure 2C). Binding and sequence 
recognition is achieved through a combination of either 



Figure 2 The DBD of human HOTl recognizes telomeric DNA in a sequence-specific manner. (A) Sequence-specific pull-down of recombinant 
HOTl, TRFl (positive control) and TBP (TATA-binding protein, negative control). Proteins were incubated with dsDNA of telomeric repeats 
(5'-TTAGGG-30, the control sequence (5'-GTGAGT-3'), the subtelomeric repeat variants (5'-TCAGGG-3', 5'-TGAGGG-3' and 5'-TTGGGG-3', as 
well as the C. elegans telomere repeat 5'-TTAGGC-3'). All DNA substrates were concatemerized from 60 bp oligonucleotides to larger DNA 
fragments (on average at least 1 kb). (B) ChIP of telomeric DNA using antibodies against HOTl, TRF2 (positive control), GFP and IgG (negative 
controls). Representative slotblot images are shown for ChIP from HeLa extracts after hybridization with a telomeric and genomic control. 
Input dilutions demonstrate the linearity of the signals acquired. (C) Structure of the DBD of HOTl bound to double-stranded telomeric DNA. 
The protein is shown as a cartoon representation (orange), whereas DNA is shown as a stick model (grey). The interacting amino acid residues 
in HOTl are shown as blue sticks, water molecules as red spheres and protein-DNA contacts are visualized as green dashed lines. 
(D) Schematic representation of all protein-DNA contacts in the complex. (E) Sequence-specific pull-down of FLAG-HOT 1 and selected 
single mutations to investigate binding specificity. Proteins were incubated with either telomeric repeats (5'-TTAGGG-3') or a control 
oligonucleotide (5'-GTGAGT-3'). (F) Atomic details of DNA sequence recognition by HOTl. K335 of helix 3 is involved in direct hydrogen 
bonding to 06 of G8 and 06 of G9. N332 of helix 3 specifically recognizes All' of the complementary strand by forming two direct H-bonds 
with the bicyclic ring system of All' (N6 and N7) (left panel). R271 of the N-terminal arm binds two bases of an AT base pair, directly to T9' and 
via a water-mediated H-bond to A12 (right panel). Source data for this figure is available on the online supplementary information page. 
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water-mediated or direct contacts with the phosphate back- 
bone (on both sides of the major groove) and the DNA 
nucleobases in the major (a3] and minor groove (N-terminal 
arm; Figures 2C and D; the structure is deposited under PBD 
entry 4J19; for crystallographic statistics please see Table I]. 



Lysine 335 is a liey residue for telomere sequence 
specificity of H0T1 

In order to bind a specific sequence, the HOTl homeo- 
domain must recognize and directly contact individual 
bases of telomeric DNA. We could identify five such 
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Figure 3 Comparison of the molecular recognition of telomeric DNA by HOTl and TRFs. (A) Schematic representation of the domain structure 
of the homeobox domains of TRFl, TRF2 and HOTl. Residues involved in DNA binding are marked with an asterisk (HOTl) or diamond 
(TRFl). Strictly conserved residues are shown with white font on red background and conserved residues are written in red font. (B) 
Superposition of structures of the HOTl DBD and TRFl DBD bound to telomeric DNA. Both binding domains recognize a different set of DNA 
bases, resulting in a different positioning relative to the 5'-TTAGGG-3' motif. 



Table I Crystallographic statistics 



Data collection 
Data set 
Beamline 
Space group 
Unit cell parameters (A) 
Wavelength (A) 
Resolution range (A)^ 
Unique reflections 
Multiplicity 
Completeness {%Y 
l/o{lY 



. (%) 



Refinement 

Resolution range (A) 

i^free (%) 
i^work (%) 

r.m.s.d. bond (A) 
r.m.s.d. angle (deg) 
B-factor 2 (A^) 

Ramachandran validation 
Favoured (%) 
Allowed (%) 
Outliers (%) 



Native 
SLS PXII 

C222i 

a= 111.4, b= 116.5, c=75.7 
1.00 
46.2-2.9 

11126 
3.6 (3.6) 
98.9 (97.4) 
18.8 (2.2) 
3.3 (45.6) 



46.2-2.90 
22.1 
17.2 
0.006 
1.0 
79.5 



99.3 
0.7 

0 



■^Values in parentheses correspond to the highest resolution shell. 

direct interactions (R271, Y327, N332, K335, R339; 
Figures 2C, D and F) in our structure. Four of these residues 
reside on helix 3, w^hich is generally referred to as the 



recognition helix and one, R271, is part of the unstructured 
N-terminal arm. 

In detail, Y327 binds to N6 of the adenine base of A6 via a 
water molecule (Figure 2F). K338 and R339 are in close 
proximity and oriented towards the phosphate of 08 and the 
nucleobase of 07, respectively, but their electron density 
was not as clearly defined, probably due to conformational 
flexibility. However, in our mutant analysis the single mutants 
K338A and R339A no longer bind DNA, and Y327A shows 
weakened binding — a clear indication that these residues are 
important for binding (Figure 2E). In addition, N332 makes 
two direct hydrogen bonds with the adenine base of All' of 
the complementary strand, but mutating N332 to A had no or 
httle effect in our binding assay. 

The interaction of a lysine, K335, of this helix might be the 
most noteworthy. The s-amino group of K335 makes a 
bifurcated hydrogen bond with carbonyl oxygens of two 
adjacent guanine bases (08 and 09) in the major groove. 
Mutation of this lysine to alanine (K335A) not only 
completely abrogates binding of HOTl to the 5'-TTAOOO-3' 
motif, suggesting that this residue is essential for binding 
(Figures 2E and F), but this mutation also results in the 
'gained' abihty to bind a nontelomeric control sequence 
(5'-OTOAOT-3') in our pull-down experiments. To our under- 
standing, this 'gain-of-function' is a strong indication that 
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K335 is a critical determinant for recognizing and binding to 
the telomeric sequence and, thus, for the specificity of HOTl. 

Cooperatively, residues in helix 3 recognize a 5'-AGGGT-3' 
( = A'] motif in the telomeric DNA. Furthermore, residues 
K325, N328 and R334 of this helix bind to the ribose 
phosphate backbone either directly (K325 to All', and R334 
to A6] or through water (N328 to All' and C12']. Apart from 
helix 3, we observe sequence-specific interactions from 
residues in the al-a2 loop (Y291 binding to the phosphate 
0 of 07] and the N-terminal arm. In full-length HOTl, this 
arm is predicted to be part of an unstructured stretch of ~ 65 
residues between the homeodomain and the preceding puta- 
tive POUs domain. Notably, in the structures of TRFl and 
TRF2 this arm was discussed to be important for extending 
the DNA-binding interface by adding specificity and affinity 
(Oourt et al, 2005). In our structure, the hydroxyl group of 
S270 of this arm binds to the 014 backbone phosphate, and 
R271 specifically contacts both base moieties of the A12-T9' 
base pair, either directly (T9'], or through a water molecule 
(A12; Figure 2F]. The importance of R271 was also confirmed 
by the fact that mutating this residue to alanine (R271A] 
completely abolishes DNA binding (Figure 2E}. As outlined in 
Figure 2D, the homeodomain of HOTl docks onto telomeric 
DNA by specifically interacting with nucleobases and back- 
bone groups of more than one telomeric repeat. In particular, 
R271 binds to the A12-T9' base pair of the following repeat. 
This overlapping binding, which was also seen in the struc- 
tures of TRFl and TRF2, accounts for the repetitive nature of 
telomeric DNA and probably provides an extra level of 
selectivity. 

While TRFl and TRF2 also recognize telomeric dsDNA via 
their respective homeobox domains, they only share low 
sequence similarity with the homeodomain of HOTl 
(13 and 5%, respectively. Figure 3 A] and DNA binding is 
different in both cases. TRFl and TRF2 both bind DNA via 
highly conserved interactions, involving direct interactions of 
a3 to TTxOO (direct contacts with phosphate (x) or bases), 
or xAOOOTx when water-mediated contacts are included. 
In HOTl, a3 of the homeodomain recognizes AxGGT directly 
(or AxGGTxAx with all resolved contacts; Court et al, 2005). 
When aligning and superimposing the DNA moieties of the 
TRFl and HOTl crystal structures, the difference in binding is 
easily visualized (Figure 3B). With respect to TRFl, HOTl is 
shifted 'down' (5' ^3'] the major groove by about one base 
and towards the subsequent telomeric repeat. 

Our crystal structure reveals in atomic detail how HOTl 
specifically recognizes and binds 5'-TTAGGG-3' repeats by 
means of its homeodomain. We could underpin these find- 
ings with mutational analyses and DNA pull-down assays, 
and we conclude that K335 is a key residue for telomere 
recognition and binding in HOTl. 

HOT 1 -telomere associations correlate with the degree 
of telomere processing 

To investigate the interaction between HOTl and telomeres 
in vivo, we analysed HOTl intracellular localization by 
immunoFISH (immunohistochemistry combined with fluores- 
cence in situ hybridization) microscopy in HeLa cells. In 
agreement with previous observations (Dejardin and 
Kingston, 2009], HOTl showed a nuclear-punctuated 
localization pattern, but in contrast to TRFl and TRF2 only 
associated with a subset of telomeres. Colocalization between 
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HOTl and telomeric DNA was observed in about 90% of all 
cells, with on average 4.5 HOTl and telomeric foci 
colocalizing per cell (Figure 4A, Supplementary Figure S4a] . 
To study the HOTl localization in another cell type, we 
performed IF stainings in mouse ES cells (Figure ID). 
Interestingly, in this setting we observed a higher degree of 
colocalization of HOTl with telomeres, with an average of 
13.5 HOTl and TRFl foci colocalizing (Figure 4B]. In telo- 
merase-positive cancer cells, such as HeLa, only a small 
fraction of telomeres is actively processed, as telomerase is 
only associated with a few telomeres at any given time (Jady 
et al, 2006; Tomlinson et al, 2006). Here most telomeres are 
believed to exist in a 'closed' state, while only some telomeres 
are in an 'open', actively processed state. In contrast, ES cells 
have high telomerase activity and maintain long telomeres 
(Varela et al, 2011). We reasoned that more colocalization 
events of HOTl with telomeres might reflect the degree of 
active telomere extension in these cells. If this was true, then 
cells that very actively extend their telomeres should have the 
highest frequency of HOTl -telomere associations. To test this 
hypothesis, we carried out IF stainings for HOTl on pachytene 
chromosome spreads from mouse spermatocytes, in which 
telomerase very actively elongates telomeres (Eisenhauer 
et al, 1997; Kozik et al, 1998; Yashima et al, 1998; Riou 
et al, 2005; Tanemura et al, 2005), using SYCP3 as a marker 
for the synaptonemal complex/the chromosome axis 
(Adelfalk et al, 2009). In agreement with a correlation 
between telomere processing and HOTl -telomere association, 
HOTl robustly localized to telomeres on all spreads analysed 
with on average 87% of chromosome ends being HOTl 
positive (Figure 4C) . We also confirmed this frequent associa- 
tion of HOTl with telomeres in mouse testes sections using 
the same IF set-up. Again, we could confirm the prominent 
localization of HOTl to telomeres in SYCP3-positive pachy- 
tene cells, highhghting the in vivo relevance of our discovery 
(Supplementary Figure S5). Together, this data vaHdates the 
in vivo localization of HOTl to telomeres and suggests that 
HOTl is a dynamic telomere-binding protein, putatively 
associated with actively processed telomeres. 

H0T1 interacts with the active telomerase complex 
and associates with CBs 

To gain further insights about HOTl function at telomeres, we 
performed HOTl immunoprecipitation assays combined with 
quantitative MS (Vermeulen et al, 2008), using both a 
polyclonal rabbit and a monoclonal mouse HOTl antibody 
and SILAC-labelled nuclear protein extracts from HeLa cells. 
HOTl itself was recovered with SILAC ratios indicative of 
specific binding to the antibody (Figure 5A and Supplemen- 
tary Figure S6a) . We also identified several proteins relevant 
to telomere homeostasis: the four core components of box 
H/ACA snoRNPs (DKCl, GARl, NHP2 and NOPIO; all part of 
the active telomerase RNP (Wang and Meier, 2004)) and the 
Ku70-Ku80 heterodimer proteins. We validated several of 
these interactions by co-IP experiments, confirming a physi- 
cal interaction of endogenous HOTl with components of 
the active telomerase complex (Figure 5B, Supplementary 
Figure S6 and Supplementary Tables S3-S6). In addition, CB 
proteins, notably the CB-scaffolding and -marker protein 
Coilin (Cioce and Lamond, 2005), were strongly enriched in 
our IP experiments (Figures 5A, B, Supplementary Figure S6a 
and Supplementary Tables S3-S6), indicating that HOTl 
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Figure 4 The degree of HOT 1 -telomere association varies between cell types. (A) Colocalization analysis of telomeres and HOTl in HeLa cells 
by immunoFISH staining. A representative image illustrating the colocalization between several HOTl foci (green) and telomeres (red) is 
shown. DAPl (blue) is used as nuclear counterstain. Colocalization events are indicated by arrows. Scale bars represent 5 |im. The 
quantification of the frequency of colocalization events was done after a 3D reconstruction of the acquired Z-stacks (^= 147). The average 
value is indicated by a red bar. (B) Colocalization analysis of TRFl and HOTl in mouse ES cells by IF staining. To visualize TRFl a LAP cell line 
(Poser et al, 2008) was used, expressing GFP-tagged TRFl at endogenous expression levels. A representative image illustrating the 
colocalization between several HOTl foci (green) and TRFl (red) is shown. DAPI (blue) is used as a nuclear counterstain. Colocalization 
events are indicated by arrows. Scale bars represent 5 |im. The quantification of the frequency of colocalization events was done after a 3D 
reconstruction of the acquired Z-stacks {_n = 126). The average value is indicated by a red bar. (C) IF stainings of HOTl at chromosome ends of 
mouse pachytene chromosome spreads. Representative images illustrating the localization of HOTl and TRF2 (in green) to chromosome ends 
are shown. The synaptonemal complex/chromosome axis is marked by SYCP3 (red) . The same field of view for the DNA counterstained by 
DAPI (greyscale) is shown in the bottom right corners. Scale bars represent 5 |im. The quantification of HOTl foci at chromosome ends was 
done after a 3D reconstruction of the acquired Z-stacks (72 = 21). The average value is indicated by a red bar. 



functions in the CB, where telomerase is assembled and 
potentially brought together with telomeres (Cioce and 
Lamond, 2005; Jady et al 2006; Tomlinson et al 2006), or 
in association with CB components in the nucleoplasm. We 
could further vahdate these interactions with analogous 
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immunoprecipitations using extracts from mouse ES cells 
(Supplementary Figure S6b and Supplementary Table S7). It 
is noteworthy that in none of our HOTl IPs any of the 
shelterin components were enriched and that in agreement 
with this HOTl was not identified in a reciprocal IP of the 
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Figure 5 HOTl associates with telomerase and CB complex components. (A) Summary of SILAC-based protein-protein interactions. 
Identification and normalized SILAC ratios are indicated for HOTl (bait) and the identified interaction partner relevant for telomere biology 
from immunoprecipitation using both a rabbit and a mouse anti-HOTl antibody. (B) Validation of the MS identifications by conventional 
immunoprecipitation. Nuclear HeLa extracts were subject to immunoprecipitation with either a polyclonal rabbit anti-HOTl or an IgG antibody, 
and were immunoblotted for DKCl, Ku70 and Coilin. HOTl IPs for the coprecipitation of DKCl and Coilin were carried out in corresponding 
LAP cell lines (Poser et al, 2008) and both proteins were detected with anti-GFP antibody. FLAG-HOT 1 was used to monitor the efficiency of the 
IP and a representative blot is shown. (C) Visualization of telomerase activity enrichment by gel electrophoresis in immunoprecipitations using 
antibodies against HOTl, DKCl (positive control), TBP, YYl, STAT3, H3K4me3 and CENP-B (negative controls), as well as TRFl and TRF2, 
using extracts from HeLa cells. All antibodies are rabbit polyclonal. A representative gel image of quantitative TRAP reaction products is 
shown. Samples were loaded on two gels and run in parallel represented by a gap between gel pictures. (D) Quantification of telomerase 
activity enrichment from the immunoprecipitations in panel C. Enrichments are normalized to immunoprecipitations using an IgG control. 
Error bars represent the s.d. of three independent experiments. Enrichments for DKCl, HOTl, TRFl and TRF2 are statistically significant with 
P<0.05 (Student's Mest). (E) Colocalization analysis of Coilin and HOTl in HeLa cells by immunofluorescence staining. A representative 
image illustrating the colocalization between several HOTl foci (green) and CBs (red; staining for Coilin) is shown. DAPI (blue) is used as 
nuclear counterstain. Colocalization events are indicated by arrows. An enhanced magnification of the boxed area is shown in the bottom right 
corners. Scale bars represent 5 |im. The quantification of the frequency of colocalization events was done after a 3D reconstruction of the 
acquired Z-stacks (n= 179). The average value is indicated by a red bar. Source data for this figure is available on the online supplementary 
information page. 



shelterin component POTl (Supplementary Figure S6c and 
Supplementary Table S8), suggesting that HOTl does not 
directly interact with the shelterin complex. 

The interaction of HOTl with the active telomerase com- 
plex components box H/ACA snoRNPs subunits (Wang and 
Meier, 2004) raises the possibility that HOTl binds to 
active telomerase. To substantiate this hypothesis, we 
performed immunoprecipitation experiments followed by 
telomerase activity measurements using the quantitative 



TRAP (telomere repeat amphfication protocol) assay for 
immunoprecipitates obtained from IPs for both HOTl and 
the positive control DKCl. Indeed, we detected telomerase 
activity in both cases (Figures 5C and 50), while five other 
nuclear DNA-binding proteins (TBP, YYl, STAT3, 
Histone3[K4me3] and CENP-B) did not show any enrichment 
of telomerase activity compared to an IgG control. Antibodies 
against TRFl and TRF2 immunoprecitated telomerase activ- 
ity, putatively dependent on the previously estabhshed link 
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between the shelterin component TPPl and telomerase (Xin 
et al, 2007]. In sum, we conclude that HOTl associates with 
the active telomerase complex. 

To corroborate the physical interaction between HOTl and 
Coilin, we performed IF stainings for HOTl and Coilin in 
nonsynchronized HeLa cells. After deconvolution and 3D 
reconstruction of the IF images the colocalization between 
both proteins was analysed. In about 85% of all cells ana- 
lysed, we observed colocalization of one to seven HOTl foci 
with CBs (Figure 5E]. Remarkably, HOTl foci preferentially 
localized to the periphery of Coilin, reminiscent of previous 
findings on the association of telomerase RNA and telomeres 
with CBs (Jady et aU 2006). This data underscores the idea 
that HOTl might be associated with actively processed 
telomeres in a telomerase-positive context. 

HOTl is a positive regulator of telomere length 

The differential degree of the HOTl -telomere colocalization, 
and the association with the active telomerase complex as 
well as the reminiscent localization to the periphery of CBs, 
suggested that HOTl might be involved in telomere main- 
tenance. To investigate this, we depleted both HOTl and 
TCABl in HeLa cells with endoribonuclease-prepared siRNA 
(esiRNA; Kittler et al, 2005; Figure 6A and Supplementary 
Figure S7]. Here TCABl, as a recently described DKCl inter- 
action partner necessary for proper trafficking of TERC 
(Venteicher et al, 2009; Zhong et al, 2011], serves as a 
reference for the extent of telomere shortening that can be 
assigned as telomerase dependent. From HOTl- and TCABl- 
depleted, as well as control-transfected cells, metaphase 
spreads were prepared 3 days post transfection and telomere 
length was determined using quantitative telomeric FISH 
(Londono-Vallejo et al, 2001]. Knockdown of both HOTl 
and TCABl resulted in significant telomere shortening 
(Figure 6B and Supplementary Figure S7] as indicated by 
reduced FISH signals. Here the entire populations of FISH 
signals shifted towards weaker signals, indicating that telo- 
meres globally shortened in the absence of HOTl. In these 
experiments the distributions after HOTl and TCABl knock- 
down were indistinguishable from each other, indicating that 
their effect on telomere homeostasis is similar. Consistent 
with this finding we also observed a significant increase in 
the appearance of signal-free chromosome ends both upon 
HOTl and TCABl knockdown (Figures 6B and C, and 
Supplementary Figure S7]. We also obtained qualitatively 
similar results by measuring telomere length after HOTl and 
TCABl knockdowns 3 days post transfection by a universal. 



single telomere-elongation length analysis (STELA; 
Figure 6D-F]. Again, depletion of both HOTl and TCABl 
led to shorter telomere length with a shortening of on average 
700 and 600 bp, respectively. Given the short nature of 
telomeres in this cell line (on average, 5kb measured by 
STELA] and the substantial contribution of variant repeats 
and subtelomeric DNA to this value (Aubert et al, 2012], our 
quantitative FISH and STELA data mutually confirm HOTl as 
a positive regulator of telomere length that quantitatively 
behaves similar to the established telomerase pathway 
member TCABl. 

Next we performed telomere length measurements after 
HOTl depletion at different time points. We reasoned that if 
HOTl is influencing telomere length in a mode dependent on 
telomerase activity, its depletion should lead to a gradual loss 
of telomeric tracts. To test this, we performed telomere length 
measurements upon HOTl depletion as a time course over 
several days. Indeed, the telomeric signal was gradually 
reduced in the absence of HOTl (Figures 7A and B]. Hence, 
the telomere length changes observed due to manipulations 
of HOTl levels are in agreement with telomerase-dependent 
changes. 

In a complementary experiment, we transiently overex- 
pressed FLAG-HOTl in HeLa cells and analysed telomere 
length 3 days after transfection. Consistent with the observed 
telomere shortening upon HOTl depletion, telomeres were 
elongated upon its overexpression as indicated by an increase 
of the entire population of FISH signals (Figures 7C and D] . 
The qFISH data suggests that the telomere length regulation 
exerted by HOTl is dependent on its telomere-binding 
capacity. To test whether HOTl -dependent telomere 
elongation indeed requires DNA binding, we used a 
HOTlAHomeobox variant. In yifro-binding studies with 
recombinant HOTlAHomeobox demonstrated that the 
deletion of this domain indeed abolished binding to telomeric 
5'-TTAGGG-3' repeats (Figure 7E] and the lack of binding 
also coincides with loss of telomeric localization in vivo 
(Supplementary Figure S4b]. To examine functional conse- 
quences of DNA-binding-deficient HOTl on telomere length 
regulation, we overexpressed FLAG-HOTl AHomeobox in 
HeLa cells. In contrast to telomere elongation upon 
expression of full-length FLAG-HOTl, the expression of 
FLAG-HOTlAHomeobox did not lead to an effect on telomere 
length in comparison to the control (Figures 7C and D]. Thus, 
the HOTl -telomere interaction is essential for telomere elon- 
gation. Collectively, these findings demonstrate that HOTl 
acts as a positive regulator of telomere length. The fact that 



Figure 6 HOTl regulates telomere length similar to the telomerase pathway member TCABl. (A] Verification of HOTl knockdown efficiency 
by western blot 48 h post transfection using the corresponding LAP cell line (Poser et al, 2008] as a reporter for protein expression. The TCABl 
knockdown was evaluated by quantitative PGR 24 h post transfection. (B] Quantification of telomere length by quantitative telomeric FISH after 
transient knockdown of HOTl and TCABl. The distributions of fluorescence intensities, in arbitrary units of fluorescence (a.u.f.], of individual 
telomeres from a total of 20 metaphases per treatment are displayed; the average intensity is indicated in red. For the gene-specific 
knockdowns, changes of average telomere signal intensity relative to the RLuc (Renilla Luciferase] control are shown (left]. Representative 
FISH images are shown for each treatment and signal-free ends are indicated by arrows (right]. Examples of individual chromosomes are 
magnified and the respective chromosomes are marked by rectangles (right]. Scale bars represent 5 iim. (C] Summary of the quantification of 
signal-free ends per metaphase after gene-specific knockdown. (D] Quantification of telomere length by universal STELA after transient 
knockdown of HOTl and TCABl. The distributions of telomere length groups in kb of individual telomeres are displayed; the average length is 
indicated in red and averages are stated with the respective s.e.m. Changes of average telomere length are shown relative to the RLuc (Renilla 
Luciferase] control. (E] Raw data of STELA reactions for the quantification of telomere length. STELA products after gel electrophoresis, 
transfer and hybridization from 12 individual reactions (lanes] per treatment are shown. Samples were loaded on two gels and were run in 
parallel, represented by a gap within HOTl esiRNA lanes. (F] Quantification of the frequency of short STELA products (< 5 kb] relative to long 
STELA products ( > 5 kb] . Five kilobases were used as a cut-off value based on this being the average telomere length as determined in the RLuc 
control. Source data for this figure is available on the online supplementary information page. 
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Figure 7 HOTl acts as a positive regulator of telomere length. (A) Verification of HOTl esiRNA knockdown efficiency using a HOT 1 -LAP cell 
line (Poser et al, 2008) as a reporter for HOTl expression by western blot after 2, 3, 4 and 6 days. Cells analysed after 4 and 6 days were 
transfected twice: After the initial transfection, cells were transfected a second time on day 3 (72 h post transfection) . All samples were run on 
the same gel, irrelevant lanes were spliced out. (B) Quantification of telomere length by quantitative telomeric FISH after transient knockdown 
of HOTl. The distributions of fluorescence intensities, in arbitrary units of fluorescence (a.u.f.), of individual telomeres from a total of 15-20 
metaphases per treatment are displayed; the average intensity is indicated in red. Changes of average telomere signal intensity are shown 
relative to the RLuc (Renilla Lucif erase) control. (C) Verification of FLAG-HOT 1 and FLAG-HOT 1 AHomeobox expression by western blot 48 h 
post transfection. (D) Results of telomere length measurements after transient overexpression of FLAG-HOT 1 and FLAG-HOT 1 AHomeobox. 
The distributions of fluorescence intensities, in a.u.f., of individual telomeres from a total of 30 metaphases per treatment are displayed; the 
average intensity is indicated in red. (E) Sequence-specific pull-down of FLAG-HOT 1 and FLAG-HOT 1 AHomeobox. Proteins were incubated 
with either telomeric repeats (5'-TTAGGG-3') or a control oligonucleotide (5'-GTGAGT-3'). Source data for this figure is available on the online 
supplementary information page. 
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Figure 8 TERT binding to chromatin is dependent on HOTl. (A) Schematic of the genetrap insertion in the Hotl gene. (B) PGR genotyping. 
Genomic DNA from wild-type, heterozygous and homozygous mice were analysed by PGR. Genotypes were confirmed using three primers 
WT-F, WT-R and MUT-R. WT-F and WT-R amplify the wild-type allele (upper band), and WT-F and MUT-R amplify the genetrap allele (lower 
band) . (C) Immunofluorescence stainings of HOTl on wild-type and Hoti^'^P^-^^'^^^'^^^'^'^/^^^^P^"^^'^^^^^^^''^ MEFs (referred to as Hotl " ^ " MEFs) . 
HOTl foci (green) are absent in Hotl~^~ MEFs. DAPI (blue) is used as a nuclear counterstain. Scale bars represent 5 |im. (D) Subcellular 
fractionation analysis of wild-type, Hotl ~ ' ~ and Tert~ ' ~ MEFs. Tubulin, PGNA and Histone3 serve as loading controls and to monitor the cell 
fractionation. Source data for this figure is available on the online supplementary information page. 



overexpression of HOTl can stimulate telomere lengthening 
further suggests that similar to TERT and TERC (Cristofari and 
Lingner, 2006), HOTl is a limiting factor, at least in this 
cellular context. 

TERT binding to chiromatin is dependent on H0T1 

The combined features of HOTl — direct binding to telomeric 
dsDNA, preferential localization to telomeres in settings of 
active processing, association with active telomerase, locali- 
zation to the periphery of CBs and a positive effect on 
telomere length— suggest a putative role of HOTl in telomer- 
ase recruitment. To test this possibility we assayed TERT 
binding to chromatin in the absence of HOTl, based on a 
previous experimental set-up that contributed to the identifi- 
cation of TPPl as a telomerase recruiter (Tejera et aU 2010). 
To avoid potential convolutions of experimental results by the 
presence of residual HOTl protein, we estabhshed a HOTl 
genetrap mouse. In Hod^^fPU-^^^^^^^^-^/f^^fP^-^^^^^^^^^^ mice, 
the genetrap construct is inserted in the first intron of the 
Hotl gene. As the start codon resides in exon 2, we expected 
that protein synthesis would be entirely abrogated (Figure 8A 
and B). To verify this, we established mouse embryonic 
fibroblasts (MEFs) of the different genotypes from littermate 
embryos and performed immunofluorescence stainings with 
our HOTl antibody. While wild-type cells showed the typical 
nuclear-punctuated pattern, HOTl signals were absent in 

^^^^Gt(pU-21T)346Card/(Gt(pU-21T)346Card ^^^^ ^p-g^^^ g^^^^ 

on this data and for the sake of simplicity, we will refer to 
this as Hotl ~ '~ genotype. We then performed the chromatin 
fractionation assay and could, as expected, detect TERT on the 
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chromatin fraction of wild-type MEFs, but not on chromatin 
fractions from TeTt~ '~ MEFs. Notably, the TERT signal was also 
absent when extracts of Hotl ~ '~ MEFs were used (Figure 8D), 
reminiscent of previous data on TPPl (Tejera et al, 2010). This 
data indicate that HOTl indeed contributes to telomerase 
recruitment in vivo. Thus, HOTl as a dynamic and differential 
telomere-binding protein contributes to telomerase-telomere 
association, mechanistically underlying its function as a positive 
regulator of telomere length. 

Discussion 

We have adapted our SILAC-based DNA-protein interaction 
screening technique to repetitive sequences as a straightfor- 
ward and robust approach to identify novel interactions. In 
addition to the established telomere-binding proteins and 
their tight interaction partner, we have identified HOTl as 
the first direct telomere-binder that positively regulates telo- 
mere length. HOTl has previously been found as a telomere- 
associated factor that localizes to some telomeres in HeLa 
cells using the PICh approach (Dejardin and Kingston, 2009). 
Interestingly, this is the only one of several recent large scale 
studies searching for telomere-associated factors, which has 
identified HOTl (Dejardin and Kingston, 2009; Giannone 
et al 2010; Nittis et al 2010; Lee et al 2011), possibly 
because also in this approach telomeric DNA was used as 
bait, rather than using the shelterin components. While by 
the PICh approach several hundreds of putative telomere- 
associated factors were identified and the PICh technique is, 
in general, addressing the global composition of (telomeric) 
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chromatin, our approach is more apt to rapidly identify the 
direct DNA binders and their tight interaction partners as 
illustrated by the identification of HOTl. Thus, our in vitro 
reconstitution assay seems particularly suited to discover 
novel direct telomere binders. 

In contrast to the other direct telomere-binding proteins, 
TRFl and TRF2, HOTl does not localize to all telomeres and 
the degree of telomere association varies between cell types. 
While in situ telomerase activity is hmited in HeLa cells to a 
simple maintenance of short telomeres, it is increased in 
mouse ES cells and is very high in mouse spermatogonia 
and spermatocytes, where telomere length is reset to the very 
long telomeres observed in mature sperm (Eisenhauer et al, 
1997; Kozik et al, 1998; Yashima et al, 1998; Riou et al, 2005; 
Tanemura et al, 2005]. The fact that the nature of the HOTl 
association with telomeres seems to be more transient in 
normally cychng telomerase-positive cancer cells compared 
to an intermediate degree of association in mouse ES cells, 
and to a prominent and robust association with telomeres in 
mouse spermatocytes, is consistent with promoting 
telomerase-dependent telomere lengthening. This raises the 
possibility that HOTl selectively interacts with actively 
processed telomeres. How exactly this interaction is 
regulated and which precise step in telomere lengthening 
HOTl promotes will have to be addressed in the future. 

While HOTl localizes to a subset of telomeres together 
with shelterin components, our HOTl immunoprecipitation 
experiments failed to establish any putative association 
between HOTl and shelterin complex members. Likewise, 
we could not retrieve HOTl after a POTl-IP, which is in 
agreement with several other studies that searched for 
shelterin-associated factors either by immunoprecipitation 
or by bimolecular fluorescence complementation (Giannone 
et al, 2010; Nittis et al, 2010; Lee et al, 2011]. Our comparison 
of the cocrystal structures of telomeric DNA with the HOTl, 
TRFl and TRF2 homeodomains indicates that while HOTl is 
shifted 'down' in 5' ^3' direction by one base towards the 
following telomeric repeat, the binding sites are largely 
overlapping. Hence, it is intriguing how these proteins 
coexist at telomeres, how they compete for binding sites 
and whether the shelterin proteins are found interspersed 
on telomeres or whether there are discrete, mutually 
exclusive patches along the telomeric tracts. This might 
further contribute to answering the question how HOTl is 
selectively restricted to a subset of telomeres. 

How does HOTl fit into the current view of telomere 
homeostasis? TCABl (WDR79/WRAP53] has been shown to 
be required for proper localization of CAB box containing 
scaRNPs to CBs, including TERC (Tycowski et al, 2009; 
Venteicher et al, 2009]. Similar to HOTl, TCABl acts as a 
positive regulator of telomere length by recruiting TERC to 
CBs (Venteicher et al, 2009] . The telomere-binding properties 
of HOTl and its functional consequences, which are 
quantitatively similar to those of TCABl, suggest that HOTl 
may act downstream of TCABl. The putative mechanism 
may involve recruiting telomerase-containing CBs to 
telomeres, and thus promoting telomerase association with 
telomeres, perhaps through additional interactions involving 
TPPl (Abreu et al, 2010; Tejera et al, 2010; Zaug et al, 2010]. 
The fact that overexpression of HOTl further stimulates 
telomere lengthening suggests that HOTl, similar to TERT 
and TERC, is also a limiting factor. Therefore, the shelterin 
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and CST complexes, and the controlled levels of TERT, TERC 
and HOTl, may all act cooperatively to define a mean 
telomere length in the cell. 

Our chromatin fractionation analysis indicates that HOTl is 
required for TERT recruitment to chromatin and, therefore, 
qualifies as a telomerase recruitment factor. Similarly, TPPl 
has been estabhshed as a telomerase recruiter likely by 
(direct or indirect] interaction with TERT via its TEL patch 
(Nandakumar et al, 2012], but based on our 
immunoprecipitation data and previous studies there is no 
indication for an interaction between HOTl and TPPl. This 
discrepancy may be explained by a step-wise telomerase 
recruitment model, in which HOTl simply acts before 
TPPl. Previous reports have shown that POTl-TPPl act as 
a processivity factor for telomerase (Wang et al, 2007; Latrick 
and Cech, 2010; Zaug et al, 2010]. Enhanced in vivo 
processivity would also result in an increased residence 
time and a higher affinity to telomeres, thus explaining the 
reduction of telomerase-telomere association in TPPl knock- 
down/knockout conditions. Given that HOTl preferentially 
localizes to telomeres in settings of high in situ telomerase 
activity, one possibility is that HOTl contributes to an 'open' 
telomere state that allows/promotes telomere elongation. In 
such a model, HOTl might contribute to rendering telomeres 
accessible and/or contributes to the delivery of telomerase to 
'open' telomeres. In this scenario, TPPl would contribute 
more prominently to maintaining telomerase at telomeres 
and, thus, to processivity stimulation. Undoubtedly, the 
mechanistic contribution of both factors and potential genetic 
interactions are exciting avenues for future research. 

Materials and methods 

Cell culture 

HeLa cells (Epitheloid carcinoma, cervix] and MEFs were cultivated 
in 4.5 g/1 glucose Dulbecco's modified Eagle's medium (DMEM] 
supplemented with 10% fetal bovine serum, lOOU/ml penicillin and 
100|ig/ml streptomycin (Gibco] at 37°C and 5% CO2. For SILAC 
labelling, HeLa cells were incubated in RPMI 1640 (-Arg, -Lys] 
medium containing 10% dialysed fetal bovine serum (Gibco] 
supplemented with 84mg/l ^^C6^^N4 L-arginine and 50mg/l 
^^C6^^N2 L-lysine (Sigma Isotec or Euriso-top] or the corresponding 
nonlabelled amino acids, respectively. Rl/E murine ES cells 
were also grown in 4.5 g/1 GlutaMAX DMEM containing sodium 
pyruvate, but supplemented with 20% FBS, 50|iM 2-mercapto- 
ethanol (Invitrogen] , LIF (MPI-CBG protein expression facility], as 
well as lOOU/ml penicillin and 100|ig/ml streptomycin (Gibco]. 
For SILAC labelling of Rl/E cells, DMEM (-Arg, -Lys] medium 
containing 20% dialysed FBS (Gibco] supplemented with 40mg/l 
^^C6^^N4 L-arginine and 80mg/l ^^C6^^N2 L-lysine (Sigma Isotec 
or Euris-top] or the corresponding nonlabelled amino acids, 
respectively, was used. Cells were collected and nuclear extracts 
were prepared as described (Butter et al, 2010]. 

Telomere pull-down 

Chemically synthesized oligonucleotides (Table II; Metabion] were 
annealed and polymerized by T4 ligase (Fermentas] and biotiny- 
lated with biotin-dATP (Invitrogen] by Klenow fragment 
(Fermentas] following the manufacturer's instructions. Twenty- 
five microgram baits were immobilized on 750 |ig paramagnetic 
streptavidin beads (Dynabeads MyOne CI, Invitrogen] and subse- 
quently incubated with either 400 |ig of SILAC-labelled nuclear 
extract or 20 |il of supernatants of crude E. colt lysates (3 |ig/ml] 
in PBB buffer (150 mM NaCl, 50 mM Tris-HCl pH 8.0, 10 mM MgCl2, 
0.5% NP-40, Complete Protease Inhibitor without EDTA (Roche]] 
for 2h at 4°C on a rotation wheel. After three washes with PBB 
buffer bead fractions were pooled, bound proteins boiled at 80°C in 
1 X LDS buffer (Invitrogen] and separated on a 4-12% gradient gel 
(Novex, Invitrogen; Butter et al, 2010]. 
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Table II Oligonucleotides used for pull-down experiments 



No Sequence motif (5' ^ 3') Primer sequence (5' 3') 



la TTAGGG for TTAGGGTTAGGGTTAGGGTTAGGGTTAGGGTTAGGGTTAGGGTTAGGGTTAGGGTTAGGG 

lb TTAGGG rev AACCCTAACCCTAACCCTAACCCTAACCCTAACCCTAACCCTAACCCTAACCCTAACCCT 

2a GTGAGT for TTGACAGTGAGTGTGAGTGTGAGTGTGAGTGTGAGTGTGAGTGTGAGTGTGAGTGTGAGT 

2b GTGAGT rev AAACTCACACTCACACTCACACTCACACTCACACTCACACTCACACTCACACTCACTGTC 

3a TTAGGG for TTAGGCTTAGGCTTAGGCTTAGGCTTAGGCTTAGGCTTAGGCTTAGGCTTAGGCTTAGGC 

3b TTAGGG rev AAGCCTAAGCCTAAGCCTAAGCCTAAGCCTAAGCCTAAGCCTAAGCCTAAGCCTAAGCCT 

4a TTGGGG for TTGGGGTTGGGGTTGGGGTTGGGGTTGGGGTTGGGGTTGGGGTTGGGGTTGGGGTTGGGG 

4b TTGGGG rev AACCCCAACCCCAACCCCAACCCCAACCCCAACCCCAACCCCAACCCCAACCCCAACCCC 

5a TCAGGG for GTCAGGGTCAGGGTCAGGGTCAGGGTCAGGGTCAGGGTCAGGGTCAGGGTCAGGGTCAGG 

5b TCAGGG rev ACCCTGACCCTGACCCTGACCCTGACCCTGACCCTGACCCTGACCCTGACCCTGACCCTG 

6a TGAGGG for GTGAGGGTGAGGGTGAGGGTGAGGGTGAGGGTGAGGGTGAGGGTGAGGGTGAGGGTGAGG 

6b TGAGGG rev ACCCTCACCCTCACCCTCACCCTCACCCTCACCCTCACCCTCACCCTCACCCTCACCCTC 



Table III Oligonucleotides used for QuickChange 



No 


Cloning/BAC engineering purpose 


Primer sequence (5' 3') 


1 
2 


HOTl homeobox domain deletion, forward primer 
HOTl homeobox domain deletion, reverse primer 


CTGCGACGAGGGAGTGCCAATATTGAAGCAGCAATCCTGG 
TGCTTCAATATTGGCACTCCCTCGTCGCAGTCGGAAAGTA 



MS data acquisition 

In-gel digestion and MS analysis was performed essentially as 
previously described by Butter et al, 2010. Peptides were desalted 
on StageTips and analysed by nanoflow liquid chromatography on 
an EASY-nLC system from Proxeon Biosystems coupled to a LTQ- 
Orbitrap XL or a Orbitrap-Velos (Thermo Electron). Peptides were 
separated on a CIS-reversed phase column (15 cm long, 75 |im inner 
diameter, packed in-house with ReproSil-Pur C18-AQ 3 jim resin 
(provided by Dr Maisch)) directly mounted on the electrospray ion 
source. We used a 107-min gradient from 2 to 60% acetonitrile in 
0.5% acetic acid at a flow of 200nl/min. The LTQ-Orbitrap XL was 
operated with a Top5 MS/MS spectra acquisition method in the 
linear ion trap per MS full scan in the Orbitrap, while for the 
Orbitrap-Velos a Top 10 acquisition method with HCD fragmentation 
was used. 



MS spectrum and data analysis 

The raw files were processed with MaxQuant (Cox and Mann, 2008) 
and searched with the Mascot search engine (Matrix Science) 
against IPI human v3.37 protein database concatenated with a 
decoy of the reversed sequences. Carbamidomethylation was set 
as fixed modification, while methionine oxidation and protein 
N-acetylation were considered as variable modifications. The 
search was performed with an initial mass tolerance of 7 p. p.m. 
mass accuracy for the precursor ion and 0.5 Da for the MS/MS 
spectra obtained with CID fragmentation and 20p.p.m. for the MS/ 
MS spectra in the HCD fragmentation mode. Search results were 
processed with MaxQuant filtered with a false-discovery rate of 
0.01. Before statistical analysis, known contaminants and reverse 
hits were removed. The protein ratios of a 'forward' experiment and 
the 'reverse' experiment were plotted in R (prerelease version 
2.8.0). Only proteins identified with at least two unique peptides 
and two quantitation events were plotted for the telomere 
pull-down (quality filter) . 

Recombinant protein expression and deletion variant 
construction 

The HOTl clone was obtained from the ORFeome collection 
(IOH40784, Invitrogen). The sequence was subcloned into SLIC- 
compatible pETM44 or pETM14 vector via SLIC cloning (Li and 
Elledge, 2007; Scholz et al, 2013) and was expressed in E. coli 
Rosetta at 18°C. E. coli extracts with overexpressed recombinant 
proteins were used for binding studies on variant DNA motives. For 
crystallization, protein was purified under nondenaturing 
conditions using Ni^^ sepharose high-performance columns (GE 
Healthcare) followed by size-exclusion chromatography on 
Superdex 75 (GE Healthcare), using 150 mM NaCl and 20 mM 
Tris-HCl (pH 7.5) as a running buffer. 



Table IV Oligonucleotides used for crystallization 





Telomeric DNA oligos 




No 


for crystallization 


Primer sequence (5' 3') 


1 


Tl; forward oligo 


CTGTTAGGGTTAGGGTTAG 


2 


T2; reverse oligo 


TCTAACCCTAACCCTAACA 



The HOTl ORFeome clone was LR recombined into a gateway 
compatible pcDNA3.1 vector with N-terminal FLAG-tag (gift of 
Christian Brandts). A FLAG-tagged HOTl homeobox deletion var- 
iant was constructed by PGR amplification of the pcDNA3.1-FLAG 
vector using primers with site-specific overhangs (Table III; 
Metabion) following the QuickChange II Site-directed mutagenesis 
kit protocol (Stratagene; Zheng et al, 2004). Point mutants were 
generated similarly using the QuickChange protocol. For expression 
in ES cells, constructs were subcloned into a pCAGGS vector. The 
constructs were sequence verified using an Applied Biosystems 
3730 Genetic Analyzer (Applied Biosystems), according to the 
manufacturer's instructions. 



Crystallization and structure determination 

The complex of human HOTl DBD and double-stranded telomeric 
DNA (5'-TTAGGGTTAGGGTTA-3') was obtained as follows: 
synthetic single-strand complementary oligonucleotides (Table IV; 
Eurofins MWG) were dissolved in water and mixed to yield a 5-mM 
stock solution. For double-strand annealing, this solution was 
heated to 96°C for lOmin and allowed to cool down to room 
temperature (RT) in the switched-off heat block within 2-3 h. This 
solution was mixed in a 1-1.2 M excess with a solution of purified 
HOTl DBD at lOmg/ml and incubated on ice for lOmin. Crystals 
of the complex were grown at 4°C by sitting-drop vapour diffusion 
from drops formed by mixing equal volumes of the complex 
solution and the crystallization condition 24% PEG 3350, 100 mM 
Na Acetate (pH 4.5) and 100 mM KSO4. For data collection, crystals 
were cryoprotected by soaking in a mother liquor supplemented 
with 20% glycerol and was flash-frozen in liquid nitrogen. 
Diffraction data were collected at the PXII beam line of the Swiss 
Light Source (SLS; Villigen, Switzerland) and processed using 
XDS (Kabsch, 1993). The structure was solved by molecular 
replacement with the programme PHASER (McCoy et al, 2007), 
using the NMR structure of the HOTl DBD as a search model 
(PDB entry 2CUF). Model building and refinement was carried out 
with Phenix (Afonine et al, 2005) and COOT (Emsley and Cowtan, 
2004), and models were validated with the programme Molprobity 
(Davis et al, 2007). 
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The alignment of the HOTl, TRFl and TRF2 homeobox domain 
sequences was done in two steps. First, a secondary structure-based 
alignment was generated using the programme Chimera 
(www.cgLucsf.edu/chimera) followed by feeding it into the web 
application Espript (espript.ibcp.fr) for a simple sequence align- 
ment scoring and colouring. 

Antibody production 

His-MBP-tagged HOTl (pETM44 vector construct) was expressed in 
E. coli Rosetta at 18°C in a fermentation tank. Cell pellet was lysed 
with Avestin and the soluble fraction subjected to affinity purifica- 
tion using Ni-sepharose. The elution fraction was concentrated with 
an Amicon Ultra 15 concentrator column and dialysed into buffer 
containing 50 mM K2PO3, 20 mM NaCl, 10% glycerol, 1 mM TCEP 
and protein inhibitors. Purified MBP-HOTl was injected into rabbits 
for immunization and the rabbits were ultimately sacrificed. The 
antibody was affinity purified using His-MBP-HOTl and MBP 
immobilized on HiTRAP desalting columns (GE Healthcare). First, 
the serum was applied to the His-MBP-HOTl column, eluted and 
then applied to the MBP column. The flow trough was quantified 
and used for subsequent experiments. 

Purified MBP-HOTl was injected into mice for immunization and 
the mice were ultimately sacrificed. Immortalized hybrid cells 
(hybridomas) were obtained by the fusion of B cells from the spleen 
of an immunized mouse with a myeloma cell line, which itself does 
not produce antibodies, using PEG and AH selection (Sigma 
Aldrich). Hybridoma clones were generated and screened using 
the Meso Scale Discovery platform (Meso Scale Diagnostics) by 
comparing affinity to His-MBP-HOTl and His-MBP-Katanin as an 
unspecific negative control. Positive clones were subcloned by 
limiting dilution and retested using the MSD platform. Based on 
the subcloned hybridoma cell lines, antibodies were purified using 
HiTRAP protein G columns (GE Healthcare) followed by acid 
elution. 

Ciiromatin immunoprecipitation 

Cells were crosslinked with 1 % formaldehyde (Thermo Scientific) 
for lOmin at 37°C, and the reaction was stopped by adding glycine 
to a final concentration of 0.125M for 5min at RT. Nuclei were 
prepared from fixed and washed cells by homogenization in cell 
lysis buffer (5mM PIPES pH 8, 85 mM KCl, 0.5% NP-40) and 
centrifugation at 1800 g for lOmin. Finally, nuclei were lysed in 
900 |il nuclei lysis buffer (50 mM Tris-HCl pH 8, 10 mM EDTA pH 8, 
1 % SDS) and lysates were sonicated for 30 min (30 s on/30 s off) in 
a Diagenode water-bath sonicator at speed 5. Following a centrifu- 
gation at 14 000 r.p.m. for 10 min, the cleared supernatants were 
snap-frozen in liquid nitrogen and stored at - 80°C. Sonication 
efficiency was routinely monitored by DNA gel electrophoresis to 
ensure that the bulk of DNA fragments was between 100 and 500 bp. 

Sonicated chromatin containing 50 |ig DNA was diluted 10 times 
in Chip dilution buffer (16.7 mM Tris-HCl pH 8, 167 mM NaCl, 
1.2 mM EDTA pH 8, 1.1 % Triton X-100, 0.01 % SDS) and precleared 
for 2 h, rotating at 4°C, with 20 \x\ blocked beads (Dynabeads Protein 
A beads (Invitrogen) incubated for 2 h with 5 mg/ml BSA) before 
the overnight incubation with 5 |ig of antibody. The following 
antibodies were used: mouse anti-HOTl (MPI-CBG Antibody 
Facility), rabbit anti-TRF2 (NBllO-57130, Novus), rabbit IgG 
(ab37415. Abeam) or mouse anti-GFP (abl218. Abeam)). The 
bound material was recovered after a 2 h incubation, rotating at 
4°C, with 30 |il blocked beads. The beads were washed for 10 min in 
each of the following wash buffers: low-salt buffer (20 mM Tris-HCl 
pH8, 150 mM NaCl, 2 mM EDTA pH 8, 1 % Triton X-100, 0.1 % SDS), 
high-salt buffer (20 mM Tris-HCl pH 8, 500 mM NaCl, 2mM EDTA 
pH 8, 1 % Triton X-100, 0.1 % SDS), LiCl buffer (10 mM Tris-HCl pH 
8, 0.25 M LiCl, 1 mM EDTA pH 8, 1 % NP-40, 1 % Na deoxycholate) 
and twice, 5 min each, in TE. ChlPed material was eluted by two 
15 min incubations at RTwith 225 |il elution buffer (0.1 M NaHCOs, 
1% SDS). Chromatin was reverse-crosslinked by adding 18|il of 
5M NaCl and incubated overnight at 65 °C, and DNA was 
submitted to RNase and proteinase K digestion and extracted by 
phenol-chloroform. 

Purified DNA recovered by ChlP was denatured in 0.2 M NaOH by 
heating to 100°C for 10 min and spotted onto a positively charged 
Biodyne B nylon membrane (Pall, VWR). Membranes were hybri- 
dized at 42°C in 6 X SSC, 0.01% SDS, 0.1% milk with 20pmol of 
DIG-labelled telomeric C-rich LNA probe (36 bp, Exiqon). Following 



hybridization washes (twice 5 min in 2 x SSC, 0.01% SDS, and 
once 2 min in 0.1 x SSC, 0.01 % SDS), the signal was revealed using 
the anti-DIG-AP antibodies (Roche) and CDP-Star (Roche) following 
the manufacturer's instructions. Images were obtained using the 
Luminescent image analyser LAS-4000 mini (GE Healthcare). 
Following the detection of a telomeric signal, membranes were 
stripped twice for 15 min each in 0.5% SDS at 60°C and twice for 
15 min each in 0.2 N NaOH, 0.1 % SDS at 37°C. Stripped membranes 
were then hybridized at 65°C in 6 x SSC, 0.01% SDS, 0.1% milk 
with P^^-random-primed labelled total genomic DNA from HeLa 
cells. Following hybridization washes (10 min in 2 x SSC, 0.01% 
SDS, 10 min in 0.5 x SSC, 0.01% SDS and 10 min in 0.1 x SSC, 
0.01 % SDS) the signal was revealed using the phosphorimager and 
an image obtained by the Typhoon scanner. 

Immunofluorescence and immunoFISH stainings on cycling 
cells 

For immunofluorescence stainings, cells were seeded 24 h before 
the treatment on either glass coverslips (0.17 mm, assorted glass, 
Thermo Scientific) or in LabTek II chambered coverglass chambers 
(Labtek). After a brief wash with Ix PBS cells were fixed in 10% 
formalin solution (Sigma Aldrich) for 10 min at RT, followed by two 
washes with 1 xPBS + 30mM glycine. Cells were then permeabi- 
lized with 1 x PBS + 0.5 % Triton X-100 for 5 min at 4°C, followed by 
again two washes with 1 x PBS + 30 mM glycine. Afterwards, cells 
were blocked in blocking solution (1 x PBS, 0.2% fish-skin gelatine 
(Sigma Aldrich)) for 15 min at RT. Primary antibodies were diluted 
in blocking solution and incubated for 1 h at RT. The following 
primary antibodies were used: mouse anti-HOTl (MPI-CBG 
Antibody Facility, 1:1000), rabbit anti-Coilin (sc-32860, Santa 
Cruz, 1:500), and goat anti-GFP (MPI-CBG Antibody Facility, 
1:2000). Cells were washed three times for 3 min each in blocking 
solution followed by a 30-min incubation at RT with secondary 
antibodies, which were diluted in blocking solution. As secondary 
antibodies, fluorescent-labelled donkey anti-rabbit-IgC or donkey 
anti-mouse-IgC antibodies with either Alexa488, Alexa555, 
Alexa594 or Alexa647 as fluorochromes (1:500, Invitrogen) were 
used. After three final washes for 3 min each in blocking solution, 
slides were briefly rinsed in distilled water and mounted using DAPI 
Prolong Gold Antifade Reagent (Invitrogen) . In the case of Labtek 
chambers, samples were incubated for 5 min with blocking solution 
containing 1 |ig/ml DAPI, followed by one wash in blocking solution. 

If FISH stainings were combined with immunofluorescence stain- 
ings (immunoFISH) , the FISH labelling was carried out directly after 
the IF protocol. Cells were post-fixed in 10% formalin solution 
(Sigma Aldrich) for 10 min at RT, followed by three washes in 70, 90 
and 100% ethanol for 5 min each, followed by a standard telomeric 
FISH protocol as previously described. 

All images were acquired with a DeltaVision Core Microscope 
(Applied Precision, Olympus 1X71 microscope) using a 100 x /1. 4 
UPlanSApo oil-immersion objective. Z-stacks (0.2 |im optical sec- 
tions) were collected and deconvolved using softWoRx (Applied 
Precision). Z-stacks were reconstructed in 3D using Imaris 
(Bitplane), and colocalization events were determined for signals 
above the background using the colocalization function. 

eslRNA synthesis 

The detailed protocol of esiRNA production has been previously 
published (Kittler et al, 2005). Briefly, optimal regions for designing 
esiRNAs were chosen using the Deqor design algorithm (Henschel 
et al, 2004) in order to fulfill two criteria: to obtain the most efficient 
silencing trigger in terms of silencing efficiency, and to get lowest 
chances to cross-silence other genes. The most favourable 
fragments were used to design gene-specific primers (Table V) 
using the Primer3 algorithm (http://frodo.wi.mit.edu/cgi-bin/ 
primer3/primer3_www.cgi). Two esiRNAs for HOTl and one each 
for TRFl and TCABl were designed and synthesized. PGR products 
for the esiRNA production were sequenced using an Applied 
Biosystems 3730 Genetic Analyzer (Applied Biosystems) according 
to the manufacturer's instructions. All positions of sequence trace 
files were confirmed by manual inspection. 

esiRNA and plasmid transfection 

For all transfections, HeLa cells were seeded in six-well plates 
(40 000 cells per well corresponding to 9.6 cm^) dish and incubated 
overnight before transfection. For esiRNA transfection, 30|il 
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Table V 


Oligonucleotides used for esiRNA production 






No 


Gene name 


ENSEMBL ID 


Forward primer (5' ^3') 


Reverse primer (5' ^ 3') 


1 
2 
3 
6 
7 


Renilla Lucifemse 
HOTl (1st esiRNA) 
HOTl (2nd esiRNA) 
TRFl 
TCABl 


ENSG00000147421 
ENSG00000147421 
ENSG00000147601 
ENSG00000141499 


GATAACTGGTCCGCAGTGGT 
ATTTCCCAAGCAGTTGTTGC 
GTTACCTCC CTGAAAGTATA 
TGGCAACTTTAAAGAAGCAGAA 
CGAATCGAGGAGCAAGAACT 


CCATTCATCCCATGATTCAA 
CCAGATCTGACAGCTTTTTGC 
GTTGGTTAT TATTTACTGGG 
AGCTTCAGTTTCCATTTCAACA 
GGGCTGAGGACATCAGAGAA 



Please note that esiRNA-template production consists of two consecutive PCRs. In the table only the transcript-specific primers for the first 
PGR round are listed. Universal tags have to be added to all primers that allow amplification with universal T7 primers in the second PGR 
round. The universal tags are 5'-TGTAAAAGGAGGGGGAGT-3' (forward) and 5'-AGGAAAGAGGTATGAGGAT-3' (reverse). The universal T7 
primers are 5'-GGTAATAGGAGTGAGTATAGGGAGAG-3' (forward) and 5'-GGTAATAGGAGTGAGTATAGGGAGAG-3' (reverse). 



Oligofectamine (Invitrogen) were diluted in 250 |il OptiMEM 
(Invitrogen) and incubated for 5min at RT. In a separate tube, 
2 |ig esiRNA were diluted in 250 |il OptiMEM. Solutions were 
combined, mixed and incubated for 20min at RT after which the 
transfection mix was evenly distributed over the dish. 

The FLAG-HOT 1 and the homeobox deletion variant construct 
were transfected using either Effectene (Qiagen) or Lipofectamine 
2000 (Invitrogen) as the transfection reagent according to the 
manufacturers' instructions. 

Telomeric quantitative FISH 

For metaphase preparation, cells were incubated for 4h with 
200 nM nocodazole in order to induce mitotic arrest. A hypotonic 
shock was achieved in 0.03 M sodium citrate at 37°G for 40min. 
Gells were fixed in an ethanol/acetic acid solution (3:1) and washed 
three times in this fixing reagent. Metaphase spreads were obtained 
by dropping suspensions of fixed cells onto clean glass slides. 

The QFISH procedure was carried out as described, using an 
Alexa488-0-0-(GGGTAA)3 or Gy3-0-0-(GGGTAA)3 PNA probe 
(Panagene) (Londono-Vallejo et al, 2001). 

Telomeric signals were quantified using the iVision software 
(Ghromaphor) . Telomere signals were segmented manually and 
average pixel intensities from every segment were quantified. For 
each metaphase, the average background intensity was determined 
and subtracted from individual telomere signals. Statistical analyses 
were done using Student's t-test. 

The number of signal-free ends per metaphase was determined by 
manual inspection of the same metaphase images that were used for 
telomere signal-intensity quantification. Statistical analyses were 
done using Student's t-test. 

Universal STELA 

One hundred picogram of Ndel/Msel (NEB) -digested DNA were 
amplified according to the previously published Universal STELA 
PGR amplification protocol (Bendix et al, 2010) with the following 
cycling conditions: 68°G for 5 min, 95°G for 2 min, 26 cycles of 95°G 
for 15 s, 58°G for 30 s and 72°G for 12 min, and 72°G for 15 min. In 
brief, 10 ng of DNA were digested with 0.4 units of Nsdel and 0.4 
units of Msel, annealed at 16°G overnight to 4.2 |iM of both 42-mer 
(5'-TGTAGGGTGA-AGAGGAGAGAAAGGGGGTGGTGGGGAGGGGGG-3') 
and 11 +2-mer (5'-TAGGGGGGTGGGG-3') panhandle-oligos with 1.3 
units of T4 ligase (NEB) and ligated to 1 nM of telorette3 (5'-TGGTG 
GGTGGATGTGGGATGGGTAAGG-3') at 35°G over night. One hundred 
picogram of ligated DNA was then used for PGR amplification with 
0.1 |iM of Adapter (5'-TGTAGG-GTGAAGAGGAGAGAA-3') and of 
Teltail (5'-TGGTGGGTGGATGTGGGATG-3') primers by the Failsafe 
Enzyme with the FailSafe buffer premix H (Epicentre). The PGR 
products were resolved on a 1.2% agarose gel, blotted onto a nylon 
membrane by capillary osmosis and hybridized at 42 °G over night 
to a DlG-labelled telomeric G-rich LNA probe (36 bp, Exiqon). The 
blots were then washed with Maleic Acid/Tween solution and DIG 
signal was revealed with anti-DIG-AP antibody (1:20 000, Roche) 
and GDP-star solution (Roche) according to the manufacturer's 
instructions. Statistical analyses were done using Student's t-test. 

Protein immunoprecipation 

Immunoprecipitations were carried out using the Dynabeads 
Protein G or A immunoprecipitation kit (Invitrogen) . Fifty microlitre 
of beads were treated with 10|ig rabbit anti-HOTl, mouse anti- 
HOTl (both MPI-GBG Antibody Facility), rabbit IgG (sc-66931. 



Santa Gruz; 2729s, Biolabs) or mouse IgG (GhromPure, Jackson 
ImmunoResearch) in PBB buffer (150 mM NaGl, 50 mM Tris-HGl pH 
8.0, 10 mM MgGl2, 0.5% NP-40, Gomplete Protease Inhibitor-EDTA 
(Roche)) and subsequently incubated with 400 |ig HeLa nuclear 
extract (SILAG-labelled if followed by MS analysis) for 2 h at 4°G on 
a rotation wheel, followed by three washes with PBB buffer. For MS, 
bead fractions were pooled, bound proteins were eluted and sepa- 
rated on a 4-12% gradient gel (Novex, Invitrogen). For co-IP 
experiments followed by western blot, bound proteins were simply 
eluted and subjected to western blot analysis. 

Western blot 

For western blot samples were boiled in Laemmli buffer (Sigma 
Aldrich) and subjected to SDS-PAGE (NuPage 4-12% Bis-Tris gels; 
Invitrogen). Gels were blotted to nitrocellulose (Protran; Schleicher 
& Schuell), blocked in 5% nonfat milk in PBST (PBS containing 
0.1 % Tween-20) for 1 h at RT and incubated over night at 4°G with 
primary antibody. The following primary antibodies were used: 
mouse anti-GFP (Roche Diagnostics, 1:4000 dilution), mouse anti- 
DMlalpha tubulin (MPI-GBG Antibody Facility, 1:50 000 dilution), 
mouse anti-Ku70 (scl7789, Santa Gruz, 1:1000 dilution), mouse 
anti-FLAG (M2, Sigma Aldrich, 1:5000) mouse anti-PGNA (sc-9847, 
Santa Gruz, 1:1000), mouse anti Histone3 (ab, 1791, Abeam, 
1:10 000) and rabbit anti-TERT (gift from Madalena Tarsounas, 
1:2000). The next day, membranes were washed three times for 
10 min each in 5 % milk PBST and were incubated for 1 h at RT with 
secondary antibody (goat anti-mouse antibody conjugated to horse- 
radish peroxidase, Bio-Rad, 1:4000, or donkey anti-mouse IRDye 
800GW or donkey anti-rabbit IRDye 800GW, both LI-GOR Odyssey, 
1:15 000). Membranes were washed three times for 10 min each in 
PBST followed by one PBS wash. Bands were visualized with 
enhanced chemiluminescence Western Blotting Detection 
Reagents (GE Healthcare) or with the LI-GOR Odyssey imaging 
system. For detection of His-tagged HOTl, TRFl and TBP, the 
Penta-His HRP Gonjugate Kit (Qiagen) was used according to the 
manufacturer's instructions. As a molecular weight standard. 
Spectra Multicolor Broad Range Protein Ladder (Fermentas), 
Seablue 2 (Invitrogen) and MagicMark XP Western Protein 
Standard (Invitrogen) were used. 

Immunoprecipitation of telomerase activity and quantitative 
TRAP assay 

For the immunoprecipitation of telomerase activity HeLa cells were 
lysed in lysis buffer (50 mM Tris-HGl (pH 8.0), 150 mM NaGl and 
1% NP-40 supplemented with Gomplete Protease Inhibitor-EDTA 
(Roche)) for 30 min on ice followed by a 30 min centrifugation step 
in a table-top centrifuge at 4°G and 20 000 g. Per IP assay, 25 |il of 
magnetic protein G beads (Invitrogen) were used. Beads were 
washed three times with 1 x PBS before use and incubated with 
5 mg/ml BSA (in 1 x PBS) for 1 h at 4°G on a rotating wheel, while 
lysates were precleared with uncoated beads for 1 h at 4°G on a 
rotating wheel. Per IP assay, 1 mg lysate was incubated with 5 ^ig 
rabbit anti-HOTl (MPI-GBG Antibody Facility), rabbit anti-DKGl 
(ab64667. Abeam), rabbit anti-TRFl (abl423. Abeam), rabbit anti- 
TRF2 (NBllO-57130, Novus), rabbit anti-TBP (sc-273, Santa Gruz), 
rabbit anti-YYl (abl2132. Abeam), rabbit anti-STAT3 (9132, Gell 
Signaling), rabbit anti-Histone3K4tri-methylated (07-473, Upstate 
Antibodies), rabbit anti-GENP-B (ab25734. Abeam) or rabbit IgG 
(sc-66931, Santa Gruz) in PBS for 2h at 4°G on a rotating wheel. 
BSA-coated beads were added, followed by a second incubation for 
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Table VI Oligonucleotides used for quantitative PCRs 



No Gene name ENSEMBL ID Forward primer (5' 3') Reverse primer (5' -> 3') 



1 GAPDH ENSG00000111640 CAGCCTCAAGATCATCAGCA TGTGGTCATGAGTCCTTCCA 

2 TCABl ENSG00000141499 AGCCAGACACCTCCTACGTG GGTTGAAGCCACAGAAGAGC 



Table VII Oligonucleotides used for cloning and BAG engineering 



No Gloning/BAG engineering purpose Primer sequence (5' 3') 



la 


HOTl BAG stitching primers to add missing 


GGATTTAGGTGGAGAAAGGAGTGGGTGGGTGTTATGGAAAGT-TAGTTGAATGAGAATGAA 




GDS and 3'-UTR, forward primer 


TAGGGAGATG 


lb 


HOTl BAG stitching primers to add 


TATAATAGAGGATTTATGATATTGTAAAGTAGTTTTAGAGAT-AGAAGGAAGGGTGTGGTG 




missing GDS and 3'-UTR, reverse primer 


GTAGATTGAA 



2 h at 4°G on a rotating wheel. Beads were then washed once with 
PBS, twice with lysis buffer and again once with PBS and finally 
recovered in 75 |il Ghaps buffer (Ghemicon). 

The quantitative TRAP assay was carried out using GoTaq 
qPGR Master Mix (Promega) and both the TS (5'-AATGGGTGGAG 
GAGAGTT-30 and AGX primer (5'-GGGGGGGTTAGGGTTAGGGTTAG 
GGTAAGG-3') at 200 nM. The reaction was run on a Mx3000p 
real-time PGR system (Stratagene) with the following protocol: 
25°G for 20min, 95°G for lOmin and 32 cycles with 95°G for 30 s, 
60°G for 30 s and 72°G for 1 min. Statistical analyses were done 
using Student's t-test. 

For verification that signals were due to the presence of the 
characteristic TRAP ladder, samples were run on a 20% precast 
TBE gel (Invitrogen) for 2h at 200 V. For size reference, the 
GeneRuler Ultra Low Range DNA Ladder (Fermentas) was used. 
The gel was stained with EtBr for visualization. 

Quantitative real-time PCR 

For quantification of TGABl mRNA levels RNA was extracted with 
the RNeasy kit (Qiagen), including DNasel digestion and from the 
eluted RNA cDNA was synthesized using the Superscript III first- 
strand synthesis kit with oligodT according to manufacturer's 
instructions. qPGR primers (Table VI) were used at 70 nM concen- 
tration together with the Absolute qPGR SYBR green mix (Abgene) 
on a Mx3000p real-time PGR system (Stratagene). Target gene 
mRNA levels were normalized against quantification of GAPDH 
mRNA levels for housekeeping. 
Statistical analyses were done using student's t-test. 

BAC TransgeneOmics 

The following BAGs (bacterial artificial chromosomes) were used in 
this study: human HOTl RP11-789B24 (Invitrogen), DKGl RPll- 
107G18, Goilin RP23-375L19 and POTl GTD-3053M7 and mouse 
TRFl RP24-402F23 (BAGPAG Resource Genter). A LAP (localization 
and affinity purification) cassette was inserted as a G-terminal 
fusion using recombineering. Isolated BAG DNA was transfected 
and selected for stable integration as described (Poser et al, 2008). 
The BAG RP11-789B24 does not cover the entire HOTl gene and was 
complemented by insertion of a cDNA fragment covering the 
missing coding and 3'-UTR sequence (Table VII). 

Immunofluorescence stainings on testes chromosome 
spreads and testes sections 

Testes were isolated from wild-type 129T2/SvEms male mice, 
between 4-6 weeks old. To prepare chromosome spreads, the 
tunica albuginea was detached and the seminiferous tubules were 
incubated in 500 |il of 1 mg/ml Gollagenase Type I (Gibco) in PBS 
for 10 min at 32 °G. The tubules were agitated slightly, and the liquid 
was removed in order to remove a proportion of interstitial cells. A 
fresh 500 |il aliquot of 1 mg/ml collagenase was added and incu- 
bated at 32°G for 30 min, disaggregating by pipetting every 10 min. 
The cells were spun at 2000r.p.m. for 5 min at 4°G. The pellet was 
resuspended in 1 ml cold PBS and filtered through a 40 |im nylon 
membrane to create a single-cell suspension. 

Fifteen microlitre of fixation and permeabilization solution (1 % 
PFA, 5 mM sodium borate pH 8.5, 0.2% Triton X-100) were added to 



each well (0 7|im) on a 10-well glass slide (StarFrost coating, 
Engelbrecht) . An ImmEdge pen (Vector Laboratories) was used to 
circle the wells and prevent leakage. The cells were diluted 1:5 into 
cold 100 mM sucrose in 1 x PBS for 2-3 min to allow hypertonic 
swelling. One to two microlitre of cell suspension were added to 
each well. The slides were incubated in a humidified chamber for 
30 min, followed by 2-3 h of drying in a laminar hood or on the 
bench at 22°G. The slides were washed briefly three times in 0.5% 
Photo-Flo (Kodak) then once in distilled water. 

To prepare testes tissue sections, whole testes were immersed in 
O.G.T Gompound (Tissue-Tek 4583) in specimen molds (Tissue-Tek 
4566 Gyromold 15 mm x 15 mm x 5 mm) and frozen at - 80°G. 
Seven-micrometer sections were cut using a Leica GM1900 and 
placed on microscope slides (StarFrost K078; 76 x 26 mm). Sections 
were fixed using 4% formaldehyde (Sigma F8775) in 1 x PBS for 
15 min at 22°G and permeabilized using 0.15% Triton X-100 for 
10 min at 22 °G and washed twice in 1 x PBS. 

To perform immunofluorescence for both chromosome spreads 
and tissue sections, blocking solution (0.2% fish gelatin (Sigma) in 
PBS - 0.1 % Tween-20) was added for 1 h at RT. Primary antibodies 
were added for 16 h at RT. The following primary antibodies were 
used: rabbit anti-HOTl (MPI-GBG Antibody Facility, 1:1000), rabbit 
anti-TRF2 (sc-9143, Santa Gruz, 1:100) and mouse anti-SYGP3 (as 
previously described (Adelfalk et al, 2009)). The cells were washed 
3 X in blocking solution for 10 min each and secondary antibodies 
(goat anti-rabbit Alexa555 or goat anti-mouse Alexa488 at 1:500, 
Invitrogen) were added for 1 h at RT. After 3 x washing in blocking 
solution for 10 min each, Vectashield (Vecta Laboratories) containing 
1 |ig/ml DAPI was added and a glass cover slip was placed on top of 
the wells and sealed closed. Slides were stored at - 20°G. 

Images were acquired with a DeltaVision Gore Microscope 
(Applied Precision, Olympus 1X71 microscope) using a 100 x /1. 4 
UPlanSApo oil-immersion objective. Z-stacks (0.2-0.5 \xm optical 
sections) were collected and deconvolved using softWoRx (Applied 
Precision). Z-stacks were reconstructed in 3D using Imaris 
(Bitplane) and colocalization events were determined for signals 
above the background using the colocalization function. For ima- 
ging of entire tubules, separate images were acquired as one panel 
and stitched together using the stitch function in softWoRx for 
visualization. 

„^f^Gt(pU-21T)346Card/(Gt(pU.21T)346Card ^ -^^ generation Of 

MEFs 

G57BL/6-GBA-Hoti^'^P^"^^'^^^^^^^''^/+ embryos were retrieved from 
the Genter for Animal Resources and Development and are based on 
the Exchangeable Gene Trap Glones system using the exchangeable 
pU21 trap vector (Araki et al, 2009) . Live animals were retrieved by 
embryo transfer and were outcrossed for several generations against 
G57BL/6 wild-type animals. MEFs were generated from littermate 
embryos at E13.5 by crossing G57BL/6-Hoti^'^P^"^^'^^^^^^^'^/ + 
animals with each other. 

Genotyping PCR 

For genotyping of Hoti^^^P^-^^^^^^^^^^^/^^^^P^-^^'^^^'^^^^^'^ mice, primers 
were designed to amplify the wild-type and mutant allele in one 
PGR reaction, both sharing the same forward primer. The PGR was 
carried out with a reaction containing 2.5 |il 10 x PGR buffer, 0.8 |il 
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Table VIII Oligonucleotides used for genotyping PCRs 



No 


Primer name 


Reverse primer (5' -> 3') 


1 


WT-F 


CCTCCTGCTCTGTCTCCTTG 


2 


WT-R 


CCCATTCAATCCCATCAAGG 


3 


MUT-R 


CCCTTCAGTCTTCCTGTCCA 



50 mM MgCl2, 0.2 |il 25 mM dNTP mix, 2|il 10|iM WT-F primer, 
each 1 |il 10 \iM WT-R and MUT-R primers, 1 |il TaqRed polymerase 
(Bioloine) and 2 |il crude genomic DNA extracts, filled up with 
HPLC-grade H2O to a total volume of 25 |il (for primer sequences 
see Table VIII). PGR conditions were as follows: initial hotstart 
denaturation at 94°G for 3 min followed by 10 cycles with 94°G for 
30 s, 62°C for 30 s with a 0.5°C touchdown decline per cycle and 
72°G for 60 s, followed by 25 cycles with 94°C for 30 s, 57°G for 30 s 
and 72°C for 60 s, finished by a final elongation at 72°C for 5 min. 
The reaction was run on a DNA Engine Thermocycler (Bio-Rad). 
After the reaction the PGR products were checked by standard gel 
electrophoresis 



Ceii fractionation 

Gell fractionation analysis for TERT binding to chromatin was 
carried out as previously described, with minor modifications 
(Tejera et al, 2010). Per sample one milhon cells were washed 
once in cold PBS and then resuspended in 200 |il ice-cold buffer 
A+ (10 mM HEPES pH 7.9, 10 mM KGl, 1.5 mM MgGl2, 0.34 M 
sucrose, 10% glycerol, Gomplete Protease Inhibitor-EDTA (Roche)). 
Samples were incubated for 5 min on ice and centrifuged for 5 min 
at 1300 g and 4°G. The supernatant was further cleaned by centri- 
fugation for 25 min at 20 000 g and 4°G, and represents the cyto- 
plasmic fraction. The pellet was washed twice in 500 |il buffer A + 
and resuspended in 100 |il ice-cold buffer B (3mM EDTA, 0.2 mM 
EGTA, Gomplete Protease Inhibitor-EDTA (Roche)). Samples were 
incubated for 30 min on ice and centrifuged for 5 min at 1700g and 
4°G. The supernatant was further cleaned by centrifugation for 
25 min at 20 000 g and 4°G, and represents the nuclear soluble 
fraction. The pellet was washed twice in 500 |il buffer B, resus- 
pended in 100 |il Laemmh buffer (Sigma Aldrich) and sonicated 
twice for 10 min each in a water-bath sonicator. This sample 
represents the chromatin fraction. The other fractions were equally 
mixed with Laemmli buffer, and all samples were boiled at 95°G for 
5 min and separated on a 4-12% gradient gel (Novex, Invitrogen). 



Supplementary data 

Supplementary data are available at The EMBO Journal Online 
(http://www.embojournal.org) . 
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